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LONG  TERM  GOAL 

The  long-term  goal  of  our  research  program  is  to  better  understand  and  quantify  the  relationships 
between  mesoscale  dynamics,  internal  waves,  turbulence  and  topographic  features  of  shallow,  tidally- 
affected  seas. 

OBJECTIVES 

The  objectives  of  the  project  during  2008  were: 

(i)  to  analyze  the  microstructure  and  velocity  structure  of  the  bottom  boundary  layer  in  the 
western  (Chinese)  sector  of  the  East  China  Sea  (ECS),  focusing  on  the  area  influenced  by 
Changjiang  (Yangtze  River)  Diluted  Water  (CDW);  continuing  investigations  on  shear 
instability  that  leads  to  the  generation  of  internal  waves  and  turbulence  in  ECS.  and 

(ii)  to  conduct  further  observations  and  analysis  on  the  genesis  of  high-frequency  non-linear 
internal  waves  packets  and  associated  turbulence  in  the  south-eastern  part  of  ECS  near  the 
ocean  shelf  break. 

APPROACH 

1 .  Analysis  of  field  data  collected  in  2006  in  ECS  during  research  cruises  of  the  Ocean 
University  of  China  (OUC). 

2.  Conduct  of  a  new  field  campaign,  led  by  the  Korean  Ocean  Research  and  Development 
Institute  (KORDI)  in  ECS,  between  the  Cheju  Island  and  the  ocean  shelf  break. 

WORK  COMPLETED 

Field  Campaigns 

The  new  research  cruise  of  KORDI  was  conducted  during  August  20-31,  2008  in  the  same  region  of 
ECS  as  the  August  2007  cruise  (Eig.  1).  A  CTD  survey  and  microstructure  measurements  were  carried 
out  using  a  SeaBird  instrument  and  TurboMAP  microstructure  profiler  at  the  hydrographic  stations 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 
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shown  in  Fig.  1.  During  this  cruise,  an  ADCP  (RD  350  kHz)  and  two  temperature-salinity  mooring 
chains  (CTR7  with  five  T/C  sensors)  were  also  deployed  in  the  survey  area. 


Fig.  1:  The  2008  cruise  of  KORDI,  August  20-31.  The  CTD  stations  were 
accompanied  by  TurboMap  microstructure  measurements. 


The  newly  obtained  data  are  now  being  processed  and  quality  assured,  and  will  be  available  for  formal 
analysis  shortly. 

RESULTS 

A,  Microstructure  in  the  Changjiang  (Yangtze  River)  Diluted  Water 

As  has  already  been  reported  (Lozovatsky,  I.D.  and  H.J.S.  Fernando,  Mixing,  internal  waves  and 
mesoscale  dynamies  in  the  East  China  Sea,  N00014-05-1-0245,  ONR  Annual  Report,  2006; 
http://www.onr.navy.mil/  soi_tech/32/reports/po_06.asp),  our  subeontractor,  OUC,  condueted  CTD, 
ADCP,  ADV,  and  microstrueture  (MSS-60)  measurements  in  ECS  at  a  station  in  CDW  (marked  by  a 
star  in  Eig.  2;  A  =  30°49’N,  cp  =  122°56’E;  the  mean  depth  Hb  =  38  m).  Two  sueeessful  periods  of 
measurements  that  lasted  25  hours  eaeh  (two  semidiurnal  tidal  cyeles)  were  completed  during 
September  3-4  and  September  7-8.  On  September  5  -  6,  a  tropical  depression  with  winds  up  to  15 
m/s  interrupted  the  ship-based  observations. 

The  barotropie  tidal  range  during  the  first  observational  period  was  1.5  m,  increasing  to  4  m  during  the 
second  one.  The  barotropie  tidal  ellipses  calculated  using  the  OTIS  software  available  through  the 
Oregon  State  University  website  (http://www.oee.orst.edu/researeh/po/researeh  /tide/  index.html)  show 
a  substantial  increase  of  tidal  flow  amplitude  during  the  seeond  period  of  observations  (see  Eig.  3). 
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This  led  to  an  increase  of  turbulence,  specifically  in  the  bottom  boundary  layer  (BBL),  by  more  than  an 
order  of  magnitude. 


Fig.  2:  The  measurement  test  site  (CDW)  of  the  2006  OUC  cruise  (September  3  -  8). 


Fig.  3:  Barotropic  tidal  ellipses  for  two  periods  of  microstructure  measurements. 

The  contour  plots  of  the  kinetic  energy  dissipation  rate  s  and  eddy  diffusivity  =  ys/N  (Fig.  4), 
where  y=  0.2  is  the  mixing  efficiency  and  the  squared  buoyancy  frequency,  show  a  turbulent  BBL 
with  a  characteristic  ^of  about  ~  10'^  W/kg  over  ~  10  meters  above  the  bottom  (mab)  during  the  near¬ 
neap  phase  of  the  barotropic  tide  (September  3-4,  Fig.  4a).  There  is  a  pronounced  phase  lag  in  ^  with 
the  distance  from  the  seafloor,  which  has  also  been  observed  in  a  rotating  tidal  flow  in  the  north¬ 
eastern  sector  of  ECS  (Lozovatsky  et  al.  2008a). 
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The  height  of  the  turbulent  layer  (where  e>  10'^  W/kg  )  inereased  twice  (~  20  mab)  and  the  dissipation 
in  the  lower  10  mab  exceeded  lO'*’  W/kg  on  September  7-8  (Fig.  4b)  closer  to  the  near-spring  tidal 
phase.  The  observed  increase  of  BBL  dissipation  is  in  line  with  the  increase  of  the  magnitude  of  tidal 
currents  .  In  fact,  according  to  the  law-of-the-wall,  s  =  Aul  jK^  at  a  specific  height  above  the 

bottom;  here  u]  =  where  and  Cd  are  the  friction  velocity  and  drag  coefficient,  respectively, 

=  1.5  in  a  rotating  tidal  flow  (Lozovatsky  et  al.  2008b)  and  /f=  0.41  is  the  Karman  constant.  The 
correlation  between  the  dissipation  rate  and  mean  velocity  near  the  bottom  was  analyzed  using  the 
ADV  data  obtained  at  the  test  site  (see  Fig.  5a). 
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Fig.  4:  The  kinetic  energy  dissipation  rate  s  and  eddy  diffusivity  Kn  during  the  25-hr  period  of 
measurement  on  September  3  -  4  (a)  and  September  7  -8  (b).  Contours  of  a  specific  potential  density  are 
overlaid  the  color  plots  to  show  the  variation  of  stratification. 


The  cubic  power  approximation,  which  is  shown  in  the  plot,  implies  that  the  mean  value  of  C ^  was 

about  1.9  X  10'\  Conversely,  the  calculation  of  the  drag  coefficient  based  on  the  quadratic  drag  law,  is 
shown  in  Fig.  5b.  The  ADV  data  obtained  in  the  2005  OUC  cruise  in  a  reversing  tidal  flow  of  ECS 
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where  used  for  this  analysis  with  A  =  1  (Lozovatsky  et  al.  2008b).  A  high  correlation  between  the 
friction  and  mean  velocities  near  the  bottom  gives  C^=1.64xl0  \  which  is  in  good  agreement  with 
the  above  estimate  obtained  for  rotating  tidal  flow  in  the  CDW  region.  Thus,  the  amplification  of 
by  a  factor  2.2  between  September  3-4  and  September  7-8  could  have  caused  the  observed  increase  of 
s  by  about  an  order  of  magnitude. 


Fig  5:  (a)  -  The  dissipation  rate  vs.  mean  velocity  at  0.3  mat  (rotating  tidal  flow  in  a  CDW  area  of  ECS, 
September  3-4  2006);  (b)  -  The  friction  velocity  vs.  mean  velocity  at  0.45  mab  (reversing  tidal  flow  near 

the  northeastern  coast  of  China  in  December  2005). 


4  2  2 

The  amplification  of  diffusivities  Kn  in  the  BBL  was  even  more  dramatic  (from  ~  10'  to  ~  10'  m  /s) 
due  to  the  combined  effects  of  an  increase  of  s  and  decrease  of  A  .  Because  the  intensification  of 
turbulence  in  BBL  and  in  the  entire  water  column  below  the  pycnocline  was  strongly  correlated  with 
the  tidal  cycle,  we  suggest  that  the  influence  of  a  mild  storm  that  passed  the  observational  area  on 
September  5-6  was  only  short  lived.  The  stratification  in  the  pycnocline  was  quickly  reestablished  by 
advection  in  the  CDW  area.  The  influence  of  tidally-induced  advection  is  demonstrated  in  Fig.  6  by 
the  depth  averaged  records  of  T{i),  S(^t),  and  aft).  For  example,  colder  and  higher-salinity  dense  water 
(At  ~  6  -  1 1  hr  and  zlt  ~  15  -  20  hr  in  Fig.  6a)  was  advected  to  the  test  site  primarily  by  the  northwest 
directed  flow  from  the  open  sea,  while  warmer  lighter  diluted  water  (At  ~  1  -  6  hr  and  zlt  ~  1 1  -  1 5  hr 
in  Fig.  6a)  was  associated  with  the  tidal  phase  when  the  flow  was  directed  from  the  northwest 
(Yangtze  River  Delta).  The  amplitudes  of  advective  variations  of  T  and  S  on  September  7-8  (Fig.  6b) 
were  as  twice  as  large  as  that  on  September  3-4  (~  1.4°C  vs.  ~  0.5°C  and  ~  0.9  psu  vs.  0.4-0. 5  psu, 
respectively)  due  to  a  longer  excursion  of  the  near-spring  tidal  current  compared  to  the  weaker  flow 
associated  with  near-neap  tide  (see  Fig.  3). 
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Fig.  6:  Time  variations  of  the  water  elevation,  Q  and  the  depth-averaged  temperature,  salinity,  and  specific 
potential  density  over  the  25-hr  period  of  September  3  -  4  (a)  and  7-8  (b). 


Intense  tidal  internal  waves  affeeted  the  small-seale  dynamics  in  the  CDW  area  during  the  second 
observational  period  (see  the  ~  10-15  meter  displacements  of  the  isopycnal  surfaces  in  the  middle  of 
the  water  column  in  Fig.  4b).  Smaller  scale  internal  waves  are  also  evident  in  this  plot,  but  they  could 
not  be  resolved  by  profiling  measurements.  An  untimely  malfunction  of  the  upward  looking  ADCP 
prevented  detailed  analysis  of  internal-wave  activity  in  the  region,  and  thus  our  research  was  focused 
more  on  turbulence  and  currents  specifically  near  the  seafloor. 

The  downward  looking  high-resolution  (dz  =  2  cm)  ADCP  allowed  to  obtain  unique  velocity  profiles 
in  the  lower  0.5  mab  with  the  lowest  level  at  0.02  mab.  During  the  25-hour  period  of  measurements 
(September  3-4),  a  very  fine  logarithmic  layer  was  observed  specifically  at  the  segments  with  an 
increased  amplitude  of  the  tidal  flow.  In  Fig.  7,  two  examples  of  U{Q  (an  hourly  averaged  data)  that 
exhibit  the  log-layers  with  a  high  confidence  fits  =  0.98  and  0.99,  respectively)  are  shown.  The 
friction  velocities  for  these  profiles  (w.  =  0.41  and  0.49  cm/s)  were  used  to  calculate  the  dissipation 
profdes  s  =  Aul Ik^  for  0.02-6  mab  with  A  =  1.5,  which  are  shown  in  Fig.  7  c,d  by  continuous 
lines  that  intersect  or  coincide  with  the  lower  ends  of  the  MSS  dissipation  profiles  (dotted  dash  lines). 
The  almost  perfect  match  between  £mss  and  ^adcp  assures  the  applicability  of  the  law-of-the-wall  with 
=  1.5  to  the  lower  --  5  m  of  the  water  column  in  ECS  (for  m.  =  0.4  -  0.5  cm/s).  When  the  flow 

magnitude  U{Q  at  <^=  0.5  mab  dropped  during  the  tidal  cycle  below  5  cm/s,  the  skin  log-layer 
approximation  failed.  The  height  of  the  mixed  BBL  here  is  controlled  by  a  characteristic  lengthscale 
A  =  u^jN ,  which  is  about  4  -  5  m  for  A~  10'^  s’'. 

An  interesting  example  of  the  dissipation  profile  is  shown  in  Fig.  7. 
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Fig:  7.  Examples  of  the  log-layer  profiles  (a:  t  =  5  hr)  and  (b:  t  =  16  hr)  and  the  corresponding  profiles  of  e 
obtained  by  the  MSS  profiling  measurements  (dotted  dash  lines)  and  the  law-of-the-wall  calculation  using 
the  ADCP  data  and  log-layer  approximations.  A  distortion  of  the  decrease  of  the  dissipation  rate  upwards 
from  the  skin  bottom  turbulent  layer  (g)  and  the  corresponding  U( Q  on  linear  scale  (e). 
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The  lowest  point  of  the  MSS  profile  at  1.5  mab  coineided  with  the  law-of-the-wall  dissipation  rate, 
but  an  internal  souree  of  turbulenee  at  ~  10  mab  (presumably  due  to  shear  instability)  “erased”  the 
expeeted  deerease  of  near-bottom  generated  turbulence  upwards  to  the  middle  of  the  water  column. 

B,  Shear  Instability  near  a  Local  Shelf  Break  in  the  Yellow  Sea 

The  study  of  flow  instabilities  and  associated  turbulence  intermittency  in  ECS  has  been  started  (Liu  et 
al.  2008)  and  will  be  continued  in  2009.  As  a  preliminary  result,  the  cumulative  probability 
distribution  of  the  gradient  Richardson  number  Ri  above  the  turbulent  BBL  using  the  2006  OUC 
measurements  near  a  local  shelf  in  the  Yellow  Sea  (Liu  et  al.  2008)  is  shown  in  Lig.  8.  The  generation 
of  turbulence  in  the  lower  part  of  water  column  was  mostly  driven  by  shear  instability  of  tidal  currents. 
The  effectiveness  of  this  process  at  different  heights  from  the  bottom  can  be  characterized  by  the 
probability  distributions  of  Ri.  Above  the  turbulent  BBL,  sporadically  the  Richardson  number  fell 
below  its  critical  value.  In  the  upper  part  of  the  logarithmic  layer  (^=  3-8  mab),  the  probability  of 
occurring  pure  Kelvin-Helmholtz  (K-H)  instabilities  was  0.32  (see  the  cumulative  distribution  function 
F{Ri)  in  Lig.  8). 


Fig.  8.  The  cumulative  distribution  of  the  Richardson  number  in  the  depth  rang  3-8  mab. 

According  to  Strang  and  Lernando  (JLM,  428,  349-386,  2001),  in  the  range  0.36  <  Ri<  1,  the 
turbulence  is  dominated  by  a  combination  of  K-H  instabilities  and  breaking  asymmetric  internal  wave. 
In  our  case,  F{Ri  <  1)  =  0.73.  Away  from  the  log-layer  (^=  9-14  mab),  the  probability  of  Ri  <  0.36 
drastically  decreases  to  0.1,  but  F{Ri  <  1),  =  0.68,  did  not  change  much.  In  the  pycnocline  {^=  15  -  20 
mab),  F{Ri  <  0.36)  for  the  occurrence  of  K-H  instabilities  increased  to  0.15  butF(R/  <  1)  decreased  to 
0.35.  Below  the  pycnocline,  the  F{Ri)  values  were  well  approximated  by  the  log-logistic  distribution 
(see  the  example  in  Lig.  8),  with  parameters  jj.  and  cr  dependent  on  the  height  from  the  seafloor. 
Statistical  analyses  of  instabilities  as  well  as  the  dissipation  rate  are  continuing  based  on  the  ADCP, 
CTD,  and  microstructure  (Turbomap)  data  obtained  in  the  summer  pycnocline  of  ECS  during  the 
recent  cruises  of  KORDI. 
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IMPACT/APPLICATION 


Our  research  program  has  been  strengthened  by  the  international  collaboration  with  Korean  and 
Chinese  oceanographers.  Two  joint  papers  with  the  OUC  colleagues  have  been  published  in  2008  in 
the  Journal  of  Continental  Shelf  Research  and  another  paper  was  accepted  for  publication  in  the 
Journal  of  Marine  Systems.  The  Pis  are  planning  to  make  a  presentation  at  the  International 
Symposium  on  Shallow  Flows  (ISSF-2008)  in  December,  Honk  Kong  and  th  eco-PI  J.-H.  Lee 
presented  two  talks  at  the  5**'  Asia  Oceania  Geosciences  Society  Meeting  in  Busan,  Korea.  A  paper  on 
a  new  BBL  mixing  parameterization  that  depends  on  the  bottom  roughness  was  published  in  2008  in 
the  Journal  of  Deep  Sea  Research.  This  research  was  supported  by  a  previous  grant  awarded  to  the  Pis 
by  ONR. 

TRANSITIONS 

None 

RELATED  PROJECTS 

The  Co-P.I.  Fernando  is  involved  in  another  ONR  funded  project  dealing  with  laboratory 
investigations  of  surface  waves  in  coastal  zone,  their  breaking  and  interactions  with  solid  objects.  This 
project  is  funded  by  the  Coastal  Geosciences  Program  of  the  ONR. 
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